Strong coupling expansion of Baxter equation in = 4 SYM 



A.V. Belitsky 



Department of Physics, Arizona State University 
Tempe, AZ 85287-1504, USA 



Abstract 

The anomalous dimensions of single-trace local Wilson operators with covariant derivatives 
in maximally supersymmetric gauge theory are believed to be generated from a deformed non- 
compact sl{2) Baxter equation. We perform a systematic expansion of this equation at strong 
coupling in the single- logarithmic limit of large conformal spin to overcome the limitation of 
the asymptotic nature of the equation. The analysis is reduced to Riemann-Hilbert problems 
for corresponding resolvents of Bethe roots in each order of the quasiclassical expansion. We 
explicitly construct the resolvents in the lowest two orders in strong coupling and find all local 
conserved charges of the underlying long-range spin chain. 



1 Introduction 



Until recently, the bulk of our knowledge about QCD or any quantum field theory has been gained 
from methods relying on perturbative expansions with respect to a small parameter. The lack of 
nonperturbative techniques overcoming limitations of perturbative considerations was a major 
obstacle in understanding strong coupling dynamics of gauge theories. The numerical framework 
of lattice gauge theory partially resolves these problems, but one still lacks reliable analytical 
methods. A concept stepping forward in recent years as a panacea for these complications is 
integrability — an idea going back to the Liouville theorem, which states that a mechanical 
system is integrable provided it has as many conserved quantities as degrees of freedom such 
that all physical quantities can be calculated exactly. 

In four-dimensional gauge theories integrable structures were first observed in studies of multi- 
reggeon compound states [H [2] . Later analyses of one- loop anomalous dimensions of the maximal- 
helicity Wilson operators in QCD [3| HI [5] revealed that their dilatation operator is mapped in 
the multi-color limit into the Hamiltonian of the noncompact XXX Heisenberg spin chain, and its 
eigenspectrum can be computed exactly by means of the algebraic Bethe ansatz [6] . Integrability 
observed in QCD anomalous dimensions at one-loop order is a generic phenomenon of four- 
dimensional Yang-Mills theories. Theories with supersymmetries obviously inherit integrability, 
although the number of integrable sectors strongly depends on the particle content of the models 
and is enhanced for theories with more supercharges eventually encompassing all operators 
in maximally supersymmetric Yang- Mills theory [H El HOl flT] . 

The gauge/string duality for the Af = 4 super- Yang- Mills theory [12] is of the weak/strong 
coupling nature and allows one to establish exact correspondence between anomalous dimensions 
of composite Wilson operators in gauge theory and energies of string excitations on the AdSs x 
background fl3[ [H]. It was demonstrated that the classical string sigma models with anti-de 
Sitter space as a factor of the target space possess an infinite set of integrals of motion and 
therefore are integrable [12]. On the gauge theory side, this result implies that the all-order 
dilatation operator for Wilson operators should be integrable. A large amount of evidence had 
been gathered to date about survival of integrability at higher orders of perturbation theory 
culminating in a conjecture of asymptotic all-order Bethe equations [161 depending on the Yang- 
Mills coupling constant g"^ = g^^Nc/ {4:7i'^) which incorporate a phase factor encoding smooth 
interpolation from weak to strong coupling. 

Recently we have formulated an alternative framework based on long-range Baxter equations 
[T7l[l8l[T9] to eigenvalue problem of all-order anomalous dimensions in maximally supersymmetric 
gauge theory. In our present note we will address, generalizing earlier considerations of Ref. |20j . 
the question of their systematic strong-coupling expansion restricting to the noncompact s/(2) 
sector and determine the spectrum of conserved charges of the magnet in the single-logarithmic 
asymptotic in the conformal spin in the lowest two orders of inverse-coupling series. One of these 
charges coincides with the eigenvalue of the dilatation operator and in the large-spin limit defines 
the cusp anomalous dimension — an observable encoding the physics of soft-gluon emission 
[211 [22] . Recently a closed integral equation was formulated for this quantity [231 [HI [21] and 
its solution at strong coupling to the lowest few orders was found numerically in Ref. [29] and 
analytically to the first [501 ED 1221 123] and eventually at an arbitrary order in Ref. [31]. These 
findings match nicely with classical [HI [33, [20], one [SHI [3Z] and two-loop calculations [55] in 
string sigma model, and with the quantum string Bethe Ansatz analysis |39j . 
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2 Asymptotic Baxter equation 



The calculation of anomalous dimensions of arbitrary twist operators is a nontrivial task even to 
one-loop order due to a large size of the mixing matrix. However, as we pointed out above the 
problem can be overcome thanks to the hidden integrability symmetry of the dilatation operator 
in A/" = 4 super Yang-Mills theory in the large- iVc limit. In the closed holomorphic sl{2) subsector 
of the theory, spanned by twist-L local Wilson operators built from scalar fields X{0) and 
light-cone covariant derivatives (schematically) 

Oiv,L(0)=tr{D^X^(0)}, (2.1) 

the one-loop mixing matrix is mapped into the Hamiltonian of the Heisenberg magnet of length 
L and spin s = 1/2 determined by the conformal spin of the field X{0). At higher orders of 
perturbation theory, the one-loop Bethe Ansatz equations [16] or Baxter equations [18] allow for 
a consistent deformation with 't Hooft coupling g and yield spectra of anomalous dimensions 
in agreement with diagrammatic calculations. They are plagued however by asymptotic nature, 
being inapplicable when the order of perturbation theory exceeds the length L of the local Wilson 
operator in question. 

In our analyses we use the method of the Baxter equation ^0] which proved to be convenient 
in analyzing various semiclassical limits including the limit of large L and N at one-loop order 

HllSl [20l H2] . The method relies on the existence of an operator Q(m) which acts on the Hilbert 
space of the chain and is diagonalized by all eigenstates of the magnet for an arbitrary complex 
spectral parameter u. For discussion of the spectrum of integrals of motion of the magnet it 
suffices to study just its eigenvalues Q{u). Thus, the long-range sl{2) Baxter equation is written 
for the polynomial 

N 

Q{u) = Hiu - uj) (2.2) 

i=i 

with zeros determined by the Bethe roots uj = Uj{g) admitting an infinite series expansion in 't 
Hooft coupling constant g, 

(x+)^el(^+("^)-^-(""))g(u + z) + (a;-)^e^(^-("'^-^+("^))Q(M - i) = t{x)Q{u) . (2.3) 

The deformation of the one-loop chain is partially encoded in the renormalized spectral parameter 
X = x[u] = ^ (^u + a/m^ — (yf^j [16], with adopted convention for x^ = x[u ± ^]. The right- 
hand side depends on the transfer matrix with the two-dimensional auxiliary space which is a 
polynomial in x with expansion coefficients determined by the conserved charges of the chaiE0. 

The most nontrivial deviations from the nearest-neighbor magnet are accommodated in the 
dressing factors a± and B in the combination 

A±(x) =or±(x) -e(x), (2.4) 

with a± responsible for the renormalization of the conformal spin of the Wilson operators (12. ip 

.■^toOp-s^)/ VZEZV ,2.5) 



^ Notice that compared to the equation in Ref. [18j we introduced a redefined transfer matrix t, t{x) — 
A+(.+)+iA_(.-)^^^^^ whose adv 
charge in the transfer matrix [17j . 



g2^+(^ )+2^ ''t{x), whose advantage compared to t{x) is the absence of ^ terms and thus corresponding 
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and providing smooth matching of the weak- and strong-couphng expansion 



-iVT^ Qm-9t)J^i s-t 

X / — I Ul- 1. (2.6) 

Jc\,,^oo] sinh^(7rK) V + gs\J \ Axx[k - gs\ 

This phase factor can be expressed as an infinite series expansion 

TV oo oo 

ze{x) = e(x+) - e(x-) =^J2J2Y1 a'^^'^^^rM [qr{x)q,{x,) - g.(x,)g,(x)] (2.7) 

j=l r=2 s=r+l 

in terms of single-excitation charges on the spin chain 



The latter define the local integrals of motion of the chain, 

2i ( 2y~^ dt , ^ 



where t/r(^) are the Chebyshev polynomials of the second kind. The anomalous dimensions of 
the operators (12.11) are related to the Hamiltonian Q2 via 

i{g) = ^^Q2{g). (2.10) 

At strong coupling, the expansion coefficients Cr^s{g) were suggested in Refs. and [21] 

at leading, first subleading and all orders, respectively. They represent analytic continuation 
of the phase (12. 6p suitable for perturbative analyses to strong coupling. The idea of using the 
strong-coupling expansion of the scattering phase (12.71) has been recently explored in Ref. [39] 
within the framework of the asymptotic Bethe Ansatz. 



3 Quasiclassical expansion 

As we pointed our earlier, due to asymptotic character of the Baxter equation it is well de- 
fined only when the interaction range does not exceed the length of the chain. Therefore, to 
study its strong-coupling expansion we have to consider a limit which does not violate this re- 
quirement. It was established in Ref. [20], that the anomalous dimensions of operators (12. ip for 
large conformal spin N and length L at strong coupling g depend on the "hidden" parameter 
^str = g'^ / L) / LF' and for ^str S> 1 the anomalous dimension becomes insensitive to the twist 
L of the operators thus circumventing limitations of the asymptotic Baxter equation. In this 
asymptotic domain, the anomalous dimension scales as ~ In N . Therefore, our goal is to develop 
the strong-coupling quasiclassical expansion in the limit g^l, N^L':^! such that ^^tr ^ 1- 



3 



To start with, let us construct a systematic expansion at strong coupling g. To this end, we 
rescale the spectral parameters with a power of the coupling constant 

u = gu, X = gx , (3.1) 

such that the Baxter equation reads 

Q{u) Q{u) 

where x^ = x[rt ± ^] and the left-hand side depends on the rescaled transfer matrix t{x) — 
x^t{x). The magnon phase 9{x) is re-expressed in terms of the local charges 

TV iV oo oo 

^(^) = XI ^(^' = ^ XI XI X ^"-''^S) [Qr{x)qs{xj) - qr{xj)qs{x)] , (3.3) 

j=\ j=l r=2 s=r+l 

with the expansion coefficients Cr^s{g) admitting the inverse-coupling expansion 

CrM = 45 + Uri + + ■■■■ (3-4) 

The Baxter polynomial in the x— variable, 

Q{u) = ll(^-x,)ll(l-^] , (3.5) 

i=i j=i ^ 

exhibits the double covering nature of the map u ^ x^ since for each value of u there are 
two corresponding values, x and l/(4a;). While the second factor is irrelevant in perturbative 
considerations, for finite coupling g these zeros do contribute on equal footing. However, as we 
will see below, their contribution vanishes from the Riemann-Hilbert problems. 
For our subsequent discussion, we introduce two resolvents 

TV TV 

^'(^) = ^X"' ^(^) = ^X7r^r-7 TT' (3.6) 



^ ' ' ' ■ -.\{x- x^ 1 - 7^ 



where r\ = {N+-^L) , which naturally emerge in the laige-g expansion of the Baxter polynomials, 

Qiu) gv 4x2-1 '^2\\grjJ Ax^-l 

rj' f±^Y fl6x^G"{x)+G"i^) 32 x' [G'{x) - G'i^)] 



^^UJ I (4.^-1)- --_,)3— +J+^(.-T 

They are not independent though and are related to each other as follows 

- 4p _ 1 + 2{2x+ 1) 2{2x-l) • ^^-^^ 
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The phases also admit the inverse-couphng expansion 

gV TV 

^{x,y) = —M%y) + -^M^,y) + -^M^,y) + oig-^) , (3.9) 
gv g V g V 

with 

Mx,y) = E E € K°H5;)g<°)(^) - g<°)(y)^s°n^)] , (3.10) 

r=2 s=r+l 



= E E - , (3.11) 

r=2 s=r+l 

t92(x,y) = E E [^r°H^)9i°Hy) - ^.°Hy)^.°^(^)] (3.12) 

r=2 s=r+l 

+ E E c£M^.°H^)^/H^)+^.'H^)^°Hy)-^.°Hy)^.^H^)-^.^Hy)^°^^ • 

r=2 s=r+l 

These are expressed in terms of the strong-couphng coefficients [23, EZl [2H] 

42 = ^^.,.-1, (3.13) 

c« = ^ ^-(-^y^' (r-i)(^-i) .3 

TT 2'-+- (s + r-2)(s-r) ' ^ ' 

45 = g + . (^ - 1)(^ - 1) , (3.15) 

and local integrals of motion 

Ux) = -qi'\x) + \q^^\x) + 0{g-') , (3.16) 

g g^ 

where 

^.°nx) = ^^, (3.17) 

U(4x2-l)3 + (4F-l)4 + (4F-l)5y' ■ ^^-^^^ 

Finally, we assume that the resolvents and the transfer matrix admit systematic expansion in 
inverse powers of the coupling 

S{x) = So{x) + -Si{x) + \s2{x) + ... , (3.19) 

g g 

G{x) = Go(x) + -Gi(x) + ^G2(x) + ... , (3.20) 

g g 

r(x) = To(x) + -Ti{x) + \t2{x) + . . . . (3.21) 

g g 

with each terms uniformly bounded. The coefficients Tk{x) are entire functions in the complex 
plane x with an essential singularity at x = 0. Substituting these expansions into the Baxter 
equation and equating the coefficients in front of powers of g~^, find a set of equations for the 
resolvents. 



5 



3.1 Leading order 

In explicit form, the leading contributions to the phase factors are 

^^^"^^ - 2{2x-l) " 2(2x + l) " (4x2-1) ' ^^-22) 

'^oW- 2(4F-1) 2(4P - 1) + (4P - 1) ' ^^-^^^ 
so that the zero-order quasiclassical form of the Baxter equation reads 

ro(x) =2cospo(x), (3.24) 

in terms of the quasimomentum 

where j3 = Lrj. From Eq. (13.241) we find that sinpo(^) is a double- valued function on the complex 
X— plane with the square-root branching points xj obeying the condition To(^j) = ±2. It becomes 
single-valued on a hyperelliptic Riemann surface defined by gluing together two copies of the 
complex X— planes along the cuts running between the branching points ^j, S = [^i, ^2] U [^3, ^4] U 
. . . . In the quasiclassical limit, these cuts accumulate the Bethe roots and correspond to regions 
of the allowed classical motion of the system in separated variables. 

From the single-valuedness of the transfer matrix, we can write down the Riemann-Hilbert 
problem 

^q(x) = vrm , (3.26) 
where here and below the principal value is defined as 

^o(x) = \ {po{x + 20) + po{x - zO)) . (3.27) 

In this study we are interested in eigenstates possessing the minimal possible energy for a given 
total spin A^. For a given total spin A^, this trajectory is realized when all cuts but two in 
S shrink into points and the Bethe roots are located on two symmetric cuts on the real axis 
[—a, —b] U [b, a] most distant from the origin. From the point of view of separated variables, this 
means that classically all but two collective degrees of freedom are frozen and the classical motion 
is confined to the two intervals. As a consequence, the complex spectral curve gets reduced to 
the genus one curve [20]. In Eq. (I3.26p . the integer m defines the position of the interval [b,a] 
inside S and m — 1 counts how many collapsed cuts are situated to the right from the interval 
[b, a] on the real axis. The minimal value of the energy is achieved for m = 1. 

Thus, to leading order of the semiclassical expansion the function Sq{x) is determined by the 
above analytical properties and prescribed asymptotic behavior at the origin and infinity 

S'q{x ^ 0) = o{x) , S'oix ^ 00) = ^ + C(x-2) , (3.28) 

for symmetric, Sq{—x) = —Sq(x), two-cut solution. It reads 

^, I f dz Wr)(z) f 1 — 4:X^ r—z -r- ab \ , . 

So{x) = -r^i> — , — ^ - + \ . 3.29 
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where the contour Ci goes around the cut [6, a] in counterclockwise direction and 



Wo{z) = 2igr])nm - -^^—^ . (3.30) 

From the asymptotic behavior of Sq{x —>■ oo) one deduces conditions on the end-points of the 
cuts 

. ^ , =0, (3.31) 

Ci 27rz ^(P-a2)(z2_52) ^ ^ 



As follows from these conditions, the branching points of the curve, ±6 and ±a, depend on the 
ratio L/N and the coupling constant g. As was demonstrated in Ref. [20], the single-logarithmic 
asymptotics emerges for the configuration when b approaches its minimal value | so that the 
inner boundaries of two cuts [— a_ — b] and [b, a] coincide with the position of poles at |x| = | and 
the outer points ±a run away to infinity: 

a^oo, b^^ + e. (3.32) 

From Eqs. (I3.3ip we find the following parametric dependence of a and b on N, L and g 

1 Amq , , , , N , , 

^ = — ^ln4a, a=^, 3.33 
y/e L 2mg 

with the parameter e related to the string parameter introduced in Ref. [20] as follows e ~ 
The analysis of the solution fl3.29p yields in the single-logarithmic regime 



s;(S) = -5^^^:5|^, (3.34) 



or in terms of the resolvent Gq, 



(3.35) 



which agrees with Ref. [39] . 

3.2 Next-to- leading order 

Collecting the terms at the order 1/(7 in the expansion of the Baxter equation, we get the equation 
for the next-to-leading correction to the resolvent 

^ 9V Po{^) -^x^Pojx) (1) 2/3x^(4x^ + 1) gr] nix) 

= — cotpo a; +^0 + TT^cotpoix)- — - . 3.36 

2 4x^ — 1 (4x"^ — 1)'* 2 smpolx) 

The Riemann-Hilbert problem for the subleading resolvent follows from this assuming that Gi{x) 
possesses the same analytic properties as the leading Go{x). Therefore, since p'oix) and smpQ(x) 
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change sign across the cuts as a consequence of the leading order Riemann-Hilbert problems 
fl3.26p . we immediately find 



grj 

"T4x2 - 1 



Pq{x + iO) cot po{x + iO) +9o\x) 



(3.37) 



where we used the fact that while S'{x) (and G{x)) is discontinuous on the cuts x G [—a, —b] U 
[b,a], 5"(^) (and G{^)) for the refiected argument is continuous. Here 



ei^\x) 



Vl{X,Z) 



4x2 - 1 vr Jc 2-Ki 4P - 1 



(3.38) 



with the contour C wrapping around the cuts [—a, —h] U [6, a] and where 

(r — l)(2m + r) 



r= 

= {Axz) 



r=2 m=0 
2 



(2m + 2r - l)(2m + 
1 

+ 



- [{2xr\2z) 



1 1 

+ 



2m+r-l (2£)2-'-(2x)2™+'-l] 

(2x + l)(2z- 1 



In- 



(4x?-l)(x-2) \{Axz-lf A{x-zfj (2x-l)(2z + l) 
The solution to the Riemann-Hilbert problem (13.371) is 

1 



3.39) 



Gi{x)=x(b p,Wi{z 



Ci 



2711 



+ 



1 I v^(P -a2)(P -62) 



^2 l-Aahz^] ^(p _a2)(p _52) 



where 



lyi(x) = 20j^^(x) - c/r7 ^^"_ ^ Po(x + iO) cotpo(^ + ^0) 



Ax 



(3.40) 



(3.41) 



in agreement with Ref. [39] • A tedious calculation yields the next-to-leading resolvent in the 
single-logarithmic asymptotics 



GAx) 



if3 



ire we 



2x(26- l)(12x2 + 1) 6Ax^(Ax^ + l)VAb^ - 1 
, r-^ m cot I X 

(4P - 1)2^x2 - 62 



{Ax 



VAb'^ - 1 
\^aP — 62 



128x3[l - 262(4x2 + 1)] _i |^Vl6s262 - 1 

cot 



(4x2 - l)Vl6x262 - 1 



2x(462 - 1) 



2Vx^ - 62 y (4x2 - l)Vx2 -62 



^ X [2(16x^ - 1) - (46^ - 1)(1 + 24x2 ^ ^g-4)j 



+ 



(4x2 - 1) - 
32x2[2x2 - 62(4x2 + 1)] (x-6)(2x + r 



(1 + ^)^(1 + 4^) 



In 



(4x2 - 1) Vx^ - fe^ (x-F6)(2x-1) 
While, the resolvent 5'^(x) can be found from it making use of the relation (13. 8p . 



(3.42) 



4 Local integrals of motion 

For the determination of the local integrals of motion Qr{g), it suffices to use reduced resolvents 
with inner end points of the cuts collided with the poles at |x| = ^. They immediately follow 
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from Eqs. fl335l) and ([332]), 



Go{x) 
GAx) 



ip 



X 



i(3 x[47rx2 + 3(4x2 + l)ln2] 



(4.1) 
(4.2) 



Since the discontinuity of the resolvent generates the distribution of Bethe roots, we conclude from 
the sub leading order in the strong-coupling expansion that the enhanced singularity at l^l = | 
implies stronger accumulation of Bethe roots around these poles. This is to be contrasted with 
weak-coupling distributions where the poles do not visible. Then the strong-couphng expansion 
of the charges 

9 

is obtained by evaluating residues of the resolvents at x = 0, 

Qi'\9) = -^f ^x~^G,{x 
g'^V J\s\<s 2vr? 

with k = 1,2. Explicitly, they read 



L r(r-l) 
L T(r) 



(4.4) 



r-2 , Ml 
+ (2r - 1 
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27r 



(4.5) 
(4.6) 



The eigenvalue of the Hamiltonian Q2 agrees with previous calculations of the cusp anomalous 
dimension at strong coupling |30l[ni[32llMllM[Il[35l|20l[Ml[37l[3g 



5 Conclusions 

In the present paper we have suggested a systematic quasiclassical expansion of all-order Baxter 
equation in the noncompact s/(2) sector of the maximally supersymmetric Yang-Mill theory at 
strong coupling. We have focused on the single-logarithmic asymptotics in the conformal spin 
corresponding to leading order contribution in the parameter e, see Eq. f l3.33p . We found the 
generating function for all local conserved charges Q,r in the first two orders of the inverse-coupling 
expansion. For r = 2, we reproduce cusp anomalous dimension at strong coupling found by other 
techniques. 

Our approach can immediately be used to compute corrections order-by-order in e— expansion. 
It can further be employed to find subsequent terms in the inverse-coupling expansion by com- 
puting higher-order terms to the resolvent. For instance, at next-to- next-to-leading order, we 
find the following Riemann-Hilbert problem 

= -i^^ry + ^0) - ^i(^)) cotpo(^)]' + ^f^(x) + 92{x) , 
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where g2{x) is 



8/3x3(1 + 16x2 + 16^4) 1 xS'o'il) 5 x^o 
^2(a;) = 



(4P - 1)5 1536 (4x2 - 1) 384 (4x2 - 1) 

1 x(21 - 296x2 + 336x4)S'^'(i) 1 x((l + 4x2)2 - 144x2)S^'(i) 



and 



384 (4x2 - 1)3 64 (4x2-1)3 

1 x(3 + 104x2 + 48x4)5^(1) 
^96 (4x2 - 1)3 ' ^^-^^ 



Its analysis deserves, however, a separate study. 

We would like to thank Gregory Korchemsky and Arkady Tseytlin for instructive commu- 
nications and very useful comments on the manuscript. This work was supported by the U.S. 
National Science Foundation under grant no. PHY-0456520. 
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